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Abstract 

The extrapolation of small-cluster exact-diagonalization calculations is used to ex- 
amine the influence of correlated hopping on valence and metal-insulator transitions in 
the one-dimensional Falicov-Kimball model. It is shown that the ground-state phase 
diagram as well as the picture of valence and metal-insulator transitions found for the 
conventional Falicov-Kimball model (without correlated hopping) are strongly changed 
when the correlated hopping term is added. The effect of correlated hopping is so strong 
that it can induce the insulator-metal transition, even in the half-filled band case. 

PACS nrs.:75.10.Lp, 71.27.+a, 71.28.+d, 71.30.+h 

In the past decade, a considerable amount of effort has been devoted to understand- 
ing of the multitude of anomalous physical properties of rare-earth and transition-metal 
compounds. Recent theoretical works based on exact numerical and analytical calcu- 
lations showed that many of these anomalous features, e.g., mixed valence phenomena, 
metal-insulator transitions, etc., can be describe very well within different versions of the 
Falicov-Kimball model (FKM) Q. For example, it was found that the spinless FKM, in the 
pressure induced case, can describe both types of intermediate- valence transitions observed 
experimentally in rare-earth compounds: a discontinuous insulator-insulator transition for 
sufficiently strong interactions § and a discontinuous insulator-metal transition for weak 
interactions |H]. In addition, at nonzero temperatures this model is able to provide the 
qualitative explanation for anomalous large values of the specific heat coefficient and for 
extremely large changes of electrical conductivity Q found in some intermediate-valence 
compounds (e.g., in SmBg). Moreover, very recently the spin-one-half version of the FKM 
has been used successfully for a description of a discontinuous intermediate-valence tran- 
sition (accompanied by a discontinuous insulator-metal transition) in SmS || as well as 



for a description of an anomalous magnetic response of the Yb-based valence-fluctuating 



These results show that the FKM could, in principle, yield the correct physics for 
describing rare-earth compounds. On the other hand it should be noted that the con- 
ventional FKM neglects all nonlocal interaction terms and thus it is questionable whether 
above mentioned results persist also in more realistic situations when nonlocal interactions 
will be turned on. An important nonlocal interaction term obvious absent in the conven- 
tional FKM is the term of correlated hopping in which the (i-electron hopping amplitudes 
Uj between neighbouring lattice sites i and j depend explicitly on the occupancy (f^fi) 
of the /-electron orbitals. To examine effects of this term on valence and metal-insulator 
transitions in the FKM we adopt here the following general form for the nearest-neighbour 
matrix elements 



which represent a much more realistic type of electron hopping than the conventional 
hopping (tij = — t if i and j are the nearest neighbours, and zero otherwise). 

Thus the spinless FKM in which the effects of correlated hopping are included can be 
written as 



" X Ujdtdj + X i' u u)-f, + /; .r/ui; «/, + uj^ftfidfdi + fth, (2) 



where /• , fi are the creation and annihilation operators for an electron in the localized 
state at lattice site i with binding energy Ef and df , di are the creation and annihilation 
operators of the itinerant spinless electrons in the d-band Wannier state at site i. 

The first term of (2) is the kinetic energy corresponding to quantum mechanical hopping 
of the itinerant <i-electrons between the nearest-neighbour sites i and j. The second term 
is just the correlated hopping term discussed above. The third term is the on-site Coulomb 
interaction between the ci-band electrons with density = \ J2i dfdi and the localized 
/-electrons with density nf = j^Y^i w ii where L is the number of lattice sites. The last 
term stands for the localized / electrons whose sharp energy level is Ef. 



compounds Q. 
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Since in this spinless version of the FKM without hybridization the /-electron occu- 
pation number f^fi of each site i commutes with the Hamiltonian (2), the /-electron 
occupation number is a good quantum number, taking only two values: Wi = 1 or 0, ac- 
cording to whether or not the site i is occupied by the localized / electron. Therefore the 
Hamiltonian (2) can be written as 

H = E h n( w ) d t d j + E f J2 w i, ( 3 ) 

<ij> i 

where hij(w) = Uj(w) + Uwi5ij and 

tij(w) = Uj + t' i:j (wi + wj). (4) 

Thus for a given /-electron configuration w = {wi, u>2 ■ ■ ■ wl} defined on a one-di- 
mensional lattice with periodic boundary conditions, the Hamiltonian (2) is the second- 
quantized version of the single-particle Hamiltonian h(w), so the investigation of the model 
(3) is reduced to the investigation of the spectrum of h for different configurations of / 
electrons. Since the d electrons do not interact among themselves, the numerical cal- 
culations precede directly in the following steps (in the next we consider only the case 
Nf + = L, which is the point of the special interest for the mixed- valence phenomena), 
(i) Having w = {u>i, W2 ■ ■ ■ wl}, U, Er and the nearest-neighbour hopping amplitudes t 
and t' fixed, (in the following t = — 1 and all energies are measured in units of t) find all 
eigenvalues of h(w). (ii) For a given Nf = J2i w i determine the ground-state energy 
E(w,U,Ef) = Y,k=i f + EfNf of a particular /-electron configuration w by filling in 
the lowest = L — Nf one-electron levels, (hi) Find the w° for which E(w,U,Ef) has 
a minimum. Repeating this procedure for different values of U,t' and Ef, one can study 
directly the the ground-state phase diagram and valence transitions (a dependence of the 
/-electron occupation number on the /-level position Ef) in the FKM with correlated 
hopping. 

The most interesting question that arises for the FKM with correlated hopping is 
whether the correlated hopping term can change the ground state phase diagram and 
the picture of valence and metal-insulator transitions found for the conventional FKM 
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(t = — 1 and t' = 0). The nature of the ground state, its energetic and structural prop- 
erties and the correlation-induced metal- insulator transitions are subjects of special inter- 
est. For the conventional FKM these problems are well understood at least at half- filling 
(Ef = 0, rif = rid = 1/2). In this case the localized /-electrons fill up one of two sub- 
lattices of the hypercubic lattice (the charge density wave state) and the corresponding 
ground state is insulating for all U > 0. Thus, for the finite interaction strength there 
is no correlation-induced metal-insulator transition in the half-filled band case. Outside 
half-filling the situation is more complicated since there exist the analytical results only for 
a restricted class of configurations 0, ||| and exact-numerical results are limited to small 
finite clusters 0,3- In spite of these restrictions both analytical and numerical results 
yield the same ground-state phase diagram which consists of only two main domains: the 
domain of the phase separation and the domain of the most homogeneous configurations. 
The domain of the phase separation is restricted to a small region of low /-electron concen- 
trations (rif < 1/4) and weak Coulomb interactions (U < 1.2). In this domain the ground 
state are incoherent (metallic) mixtures |3|, [?j of the empty configuration (whose length I 
is at least L/2) and a configuration w (whose length is I < L/2). Outside this domain 
the ground states are the most homogeneous configurations, which are insulating for all 
/-electron concentrations. The boundary between these two domains is the boundary of 
correlation- induced metal- insulator transitions. 

One can expect, on the base of simple arguments, that the ground-state phase diagram 
of the FKM with correlated hopping will be fully different from one discussed above for the 
conventional FKM. Indeed, the following selection of hopping matrix amplitudes t = — 1 
and t' > may favor the segregated configuration {11 . . . 100 . . . 0} since the itinerant 
d electrons have the lower kinetic energy in this state. This mechanism could lead, for 
example, to the instability of the CDW state {10 . . . 10} that is the ground state for t' = 
and thereby to a metal-insulator transition, even in the half-filled band case. To verify this 
conjecture we have performed an exhaustive study of the model on finite clusters (up to 
32 sites) for a wide range of parameters t' and U. The results of numerical calculations 
are summarized in Fig. 1 in the form of the t'-U phase diagram. In addition to the 
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alternating configuration w\ = {10 . . . 10} that is the ground state at t' = for all nonzero 
U we found two new phases that can be the ground states of the model, and namely, 
the alternating configuration W2 = {1100 . . . 1100} with double period and the segregated 
configuration W3 = {11 . . . 100 . . . 0}. Thus at nonzero t' the configuration w\ becomes 
unstable against the transition to W2 and 1^3 . The transition from w% to W2 is the insulator- 
insulator transition since both w\ and wi configurations have the finite gaps A at the Fermi 
energy || for nonzero values of U (see the inset in Fig. 1). On the other hand the segregated 
configuration W3 is metallic (A = 0) for all nonzero U and so the transition from w\ to 
u>3 (as well as from W2 to W3) is the insulator- metal transition. Thus we arrive to the 
very important conclusion, and namely, that the correlated hopping term can induce the 
insulator-metal transition, even in the half-filled band case. 

The fact that the correlated hopping plays the important role in stabilizing the ground 
state indicates that the picture of valence and metal-insulator transitions found in our 
previous papers within the conventional FKM g, ^] could be dramatically changed if finite 
values of t' will be considered. To verify this conjecture we have performed the study of the 
model for two representative values of U {U = 1.5 and U = 4) and a selected set of i' values. 
The results of numerical calculations obtained for U = 1.5 on finite clusters (up to 32 sites) 
can be summarized as follows: (i) For sufficiently small values of t' (t' < 0.3) the ground 
states are the most homogeneous configurations for all /-electron concentrations, (ii) For 
intermediate values of t' (t' ~ 0.6) the ground states are the segregated configurations for 
rij < 0.5 and the most homogeneous configurations for rif > 0.5. (iii) For large values of 
t' {if > 0.8) the ground states are only the segregated configurations. Since no significant 
finite-size effects have been observed on clusters up to 32 sites we suppose that these 
results can be satisfactory extrapolated on large systems and used for a study of effects of 
correlated hopping on valence and metal-insulator transitions. Typical valence transitions 
obtained on the extrapolated set of configurations are presented in Fig. 2. For t' = 0.3 the 
valence transition holds all characteristic features found in our previous papers for t' = 0. 
The basic structure of the transition is formed by the most homogeneous configurations 
with the smallest periods (n/ = 1/5, 1/4, 1/3, 1/2, 2/3, 3/4, 4/5). This primary structure is 
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independent of L while the secondary structure corresponding to remaining configurations 
strongly depends on L and forms the gradual transition between two valence states from 
the primary structure. The correlated hopping term reduces considerably the width of 
the primary structure, especially for nj > 0.5, but it does not change the nature of the 
ground state. The ground states for t' = 0.3 are the most homogeneous configurations 
which are insulating (see the inset in Fig. 2) and thus there is no insulator-metal transition 
induced by external pressure for small values of t'. (It is generally supposed that in the 
pressure induced case Ef is proportional to pressure.) From this point of view the situation 
looks more hopeful for intermediate values of t'. As was mentioned above for this case the 
ground states are the segregated (metallic) configurations (for n/ < 0.5) and thus one 
could expect the pressure induced insulator-metal transition at Ef = E c (corresponding 
to rif = n c = 0.5), of course if the segregated configurations persist as the ground states 
on some finite region of Ef values. Fig. 2 illustrates that such a situation indeed becomes. 
For t' = 0.6 the segregated configurations are the ground states in the wide region above 
E c = —0.02 and the valence transition in this region is continuous. Below E c the ground 
states are the most homogeneous (insulating) configurations. In this region the primary 
structure is practically fully suppressed and only the intermediate valence states with 
rif = 1/2 and rif = 2/3 have the finite regions of stability. At Ef = E c the FKM exhibits a 
discontinuous insulator-metal transition that is apparently induced by correlated hopping. 
This result shows how crucial role plays the correlated hopping in the mechanism of valence 
and insulator- metal transitions. Already relatively small values of t' can fully destroy the 
picture of valence and metal-insulator transitions found for the conventional FKM and 
therefore the correlated hopping term cannot be neglected in the correct description of 
valence and metal-insulator transitions. Of course, for large values of t' we can expect 
the largest deviations from t' = case, since the ground states in this region are only 
the segregated configurations which contrary to the most homogeneous configurations are 
metallic. The inset in Fig. 2 confirms this conjecture. The valence transitions for large 
values of if are only continuous and fully differ from if = case. 

In the strong coupling region (U = 4) the ground states are either the most homoge- 
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neous configurations (for small and intermediate t') or the segregated configurations (for 
large t'). Thus for small and intermediate values of t' the strong coupling picture of valence 
and metal-insulator transitions will be the same as one for t' = 0. All ground states are in- 
sulating and the valence transitions, for both conventional FKM and FKM with correlated 
hopping consist of a few discontinuous valence transitions which number is reduced with 
increasing U. Only one difference between valence transitions obtained within these two 
models is that the width of the valence transition in the FKM with correlated hopping is 
reduced due to finite t'. Fig. 3 illustrates that the effect is very strong even for relatively 
small values of t'. This indicates that the correlated hopping plays important role also in 
the strong coupling limit although unlike the weak-coupling limit now small values of t' do 
not change qualitatively the picture of valence and metal-insulator transitions. In order 
to change qualitatively the picture of valence and metal-insulator transitions in the strong 
coupling limit, larger values of t' should be considered. For example, we have found that 
for U = 4 and t' = 1.2 the ground states are the segregated (metallic) configurations for 
all /-electron concentrations and the corresponding valence transition is continuous. 

The exactly opposite effect on valence transitions in the strong coupling limit has the 
correlated hopping term with negative values of t'. As shown in Fig. 4. the correlated 
hopping with t' < stabilizes the insulating intermediate-valence states corresponding to 
the most homogeneous configurations with the smallest periods (which are the ground 
states in this region) and this effect is very strong even for relatively small values of t' . 
This again leads to conclusion that the correlated hopping must be taken into account in 
the correct description of valence and metal- insulator transitions. 

This work was supported by the Slovak Grant Agency VEGA under grant No. 2/4177/97. 
Numerical results were obtained using computational resources of the Computing Centre 
of the Slovak Academy of Sciences. 
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Figure Captions 

Fig. 1. t'-U phase diagram of the FKM with correlated hopping at half-filling (Ef = 
0,n/ = rid = 0.5). Three different phases correspond to the alternating configuration 
w\ = {10... 10}, the alternating configuration w<i = {1100 ... 1100} with double period 
and the segregated configuration u>3 = {1 . . . 10 . . . 0}. For L = 24 the phase diagram 
has been obtained over the full set of /-electron configurations, while for L = 600 only 
the restricted set of configurations 101,102 and w% has been used to determine the phase 
boundaries. The inset shows t'-dependence of energy gaps corresponding to the ground 
state configurations for different values of U. 

Fig. 2. Dependence of the /-electron occupation number n/ on the /-level position Ef 
for different values of correlated hopping t' at U = 1.5. The upper inset shows energy gaps 
corresponding to the most homogeneous configurations that are ground states for if = 0.3 
(for all /-electron concentrations) and if = 0.6 (for n/ > 0.5). The lower inset shows 
examples of rif(Ef) behaviour for large values of if. 

Fig. 3. Dependence of the /-electron occupation number nj on the /-level position 
Ef for different values of correlated hopping if at U = 4. The inset shows energy gaps 
corresponding to the most homogeneous configurations that are ground states for all /- 
electron concentrations for if = 0.4 and U = 4 

Fig. 4. Dependence of the /-electron occupation number n/ on the /-level position 
Ef for different values of correlated hopping t' < at U = 4. The inset shows energy 
gaps corresponding to the most homogeneous configurations that are ground states for all 
/-electron concentrations for t' = —0.2 and if = —0.4. 
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